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Abstract 

The beta delayed deuteron emission from 6 He is studied in a dynamical mi- 
croscopic cluster model. This model gives a reasonably good description for all 
the subsystems of 6 He and 6 Li in a coherent way without any free parameter. 
The beta decay transition probability to the 6 Li ground state is underesti- 
mated by a few percents. The theoretical beta delayed deuteron spectrum 
is close to experiment but it is also underestimated by about a factor 1.7. 
We argue that, in spite of their different magnitudes, both underestimations 
might have a common origin. The model confirms that the neutron halo part 
of the 6 He wave function plays a crucial role in quenching the beta decay 
toward the a + d channel. 

PACS numbers: 23.40.Hc, 21.60.Gx, 27.20.+n 
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I. INTRODUCTION 



In the past few years the large amount of experimental data which have been accumulated 
on the structure and reactions of unstable nuclei revealed the existence of a neutron halo 
structure in nuclei with a large neutron excess |],|2|]. The best known of these nuclei is 11 Li. 
However, there are ambiguities in the theoretical description of this nucleus mainly because 
few 9 Li + n scattering data are available, and because the 9 Li core is soft @. Fortunately, 
these problems do not occur in the case of the other prominent neutron halo nucleus, 6 He. 
Sophisticated calculations exist which can account for several properties of 6 He |3|,(|,f)]||[7|] . 

One of the most surprising feature of 6 He is the very small Gamow- Teller (GT) beta decay 
branching ratio toward the a + d continuum channel. The (3 delayed deuteron emission from 
6 He was first observed in Ref. ||, and the measured branching ratio (2.5 ± 0.5) • 10~ 6 is two 
orders of magnitude smaller than the result of a phenomenological analysis of the decay 
process in terms of an R- matrix formalism Later an experiment with higher statistics 
resulted in a new value (7.6 ± 0.6) ■ 10~ 6 H. Even this larger value is overestimated in a 



potential model JXO] and in a three-body a + n + n model [11 



A semimicroscopic model [0 shows that the smallness of the branching ratio is the 
result of a cancellation, taking place between two parts of the GT matrix element, which 
have different signs. The 'internal' part comes from the typical nuclear regions (r < 5 fm) 
of the 6 He and a + d wave functions. The 'external' part comes from the 5 fm < r < 20 
fm region. The importance of the latter part is a clear consequence of the halo structure 
in 6 He. In that semimicroscopic model, the 6 He initial state is described with a 6-nucleon 
microscopic wave function but the final state is treated in a potential model. This hybrid 
treatment is flexible enough to reveal the quenching mechanism of the branching ratio. 
However, its agreement with experiment is probably partly fortuitous because the results 
are very sensitive to model assumptions. In Ref. jnj the branching ratio is overestimated 
by an order of magnitude because the assumption of a pure a + dineutron configuration 



enhances the 'external' part of the GT matrix element. In Ref. |TT[ the overestimation is 
reduced because of a much more realistic description of 6 He but the treatment of the a + d 
scattering state is not consistent with the treatment of the ground state of 6 He. 

Clearly, an accurate and consistent description of this process is desirable. However, 
obtaining reliable results requires that the following conditions are met. The model must 
reproduce in a realistic way (i) the deuteron binding energy and size, (ii) the low-energy 
n + n phase shifts, (iii) the a + n scattering, (iv) the a + d scattering, and (v) the 6 He 



binding energy and size. In addition the analysis of Ref. [Iz] indicates that (vi) the halo 
must be well described up to distances as large as 20 fm. In this paper we calculate the j3 
delayed deuteron spectrum in a fully microscopic model. We try to use a parameter free 
model without any ad hoc model assumption. With the Minnesota interaction conditions 
(i) to (iv) are fulfilled [0,0] and a good 6 He wave function satisfying (v) is available J7| 
although the basis has to be enlarged to meet (vi). The fact that a few percent error in the 
model at typical nuclear scales can modify the order of magnitude of the final result, as is 
pointed out in Ref. | I2|| , calls for a test example which is sensitive to the details of the wave 
functions only up to less than 10 fm. Our test case will be the (3 decay transition to the 
ground state of 6 Li. 

The model is recalled in Sec. II. The results are described and discussed in Sec. III. 
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Conclusions are presented in Sec. IV. 



II. MODEL 

The (3 delayed deuteron emission probability per time and energy units, dW/dE, can be 
expressed [|T^,|I^] as 

% = ^G*f(Q-E)B GT (E), (1) 
dE ir^vn p 

where m is the electron mass, v is the relative velocity between a and deuteron and Gp = 
3.002 -i0~ 12 is the dimensionless /3-decay constant. The phase space factor, or Fermi integral, 
/ depends on the kinetic energy Q — E available for the electron and antineutrino. The mass 
difference Q between the initial and final particles is 2.03 MeV. The maximum total energy 
(including mass energy) is Wo = 2.54 MeV. The GT reduced transition probability reads 

Bgt(E) = £ | (W« M , | £ ^ | *S) | 2 , (2) 

\ ZJ "I" L ) MM'n j=l 

where ty ad and \l/ 6He are the microscopic wave functions of the a + d scattering state and of 
the 6 He ground state, respectively (we shall specify them later), P and cr J are respectively 
the isospin operator and Pauli spin-matrices of nucleon j, and A = 1.26 is the ratio of the 
axial- vector to vector coupling constants. The half life tyz f° r the /3-decay to the ground 
state of 6 Li is calculated with 

(/V)- 1 = (21n2)7r 3 (mc 2 A)G 2 J B GT (g.s.), (3) 

where Bcr(g.s.) is given by 

£ GT (g.s.) = —f— £ | Y.±ol I $S)| 2 . (4) 

In a microscopic dynamical model has been developed for the description of the ground 
state of 6 He. In the present work we use the wave function of that model 



*S = E E + 

S,h,l2,L S,l\,l2,L 

Na-1 

= E E ^{[^n^l} S ^SL(Pnn,Pa { nn))} JM } 

S,h,h,L i=0 

N a -1 

+ E E ^{[r(^^kS(p Qfll p„ w )] JM } 
+ T,^{[i^%x t Kpu)} m }- (5) 



Here ^4 is the intercluster antisymmetrizer, the p vectors are the different intercluster Jacobi 
coordinates, and [ ] denotes angular momentum coupling. While $ n is a neutron spin-isospin 
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eigenstate, $* is the antisymmetrized triton internal state in the harmonic oscillator shell 
model with a single oscillator parameter. The antisymmetrized ground state (i = 0) and 
monopole excited states (i > 0) of the a particle are represented by the wave functions 



N a 



*? = 2J^ft> i = Q,l,...,(N a -l), (6) 

where <f>p is a translation invariant shell-model wave function of the a particle with size 
parameter j3j and the A{j parameters are to be determined by minimizing the energy of the 
a particle |16|]. Similarly the size parameter of the tritons is determined from the energy 



stabilization. We choose the same parameters for the wave function as in Ref. J7| with N a 
3. The following SL terms are included in the S, {lil-^L coupling scheme for the J n = + 
ground state of 6 He: {a(nn); 0(00)0}, {a(nn); 1(11)1}, {n(an); 0(00)0}, {n(an); 0(11)0}, 
{n(an); 1(11)1}, and {tt,00}. Putting (|5|) into the six-nucleon Schrodinger equation, we 
arrive at an equation for the intercluster relative motion functions X- These functions are 
expanded in terms of products of tempered Gaussian functions exp(— 7»p 2 ) [T7|] with different 
ranges ji for each type of relative coordinate. The expansion coefficients are determined from 
a variational principle. 

The Minnesota force |13j with an exchange parameter u = 0.92 and a slightly modified 
spin-orbit component (see Ref. 0) reproduces very well the experimental phase shifts of 
all the N + N and a + N scattering states which take place in (|5|) 0. After this force 
choice there remains no free parameter in the model. This model gives for 6 He an a + n 
+ n three-cluster separation energy of 0.961 MeV, i.e. practically the experimental value 
0.975 MeV. It was pointed out in Ref. @ that this value could not be reproduced without 
the t + t component. As the GT matrix element is sensitive to the details of the 6 He wave 
function up to 15-20 fm in the p a ( nn ) coordinate [[L2[], our wave function should be correct 
in this region. To achieve this, we enlarge the number of basis states which describe large 
Pa(nn) separations in such a way that the range of our last basis functions is ~ 20 fm in these 
relative motions. 

To be consistent, we must employ the same model parameters in the description of the 
a + d bound and scattering states as for 6 He. For example the use of a different nucleon- 
nucleon interaction for a + d would affect the product of the spatial parts of the 6 He and 
a + d wave functions and modify the balance between the different parts in the GT matrix 
element, and could lead to inaccurate results. The Minnesota force reproduces the N + N 
effective range parameters in the triplet even partial wave by assuming a pure 3 Si state, i.e. 
without tensor coupling with a D state. However, if a tensor force were taken into account 
in addition to the Minnesota force, it would play a very large role in this partial wave in 



a coupled S\- D x description |L8|. This means that if a pure 6 Si state reproduces the 



experiment, then the effective interaction in this partial wave is too strong. As was pointed 
out in 0, the Minnesota interaction overbinds the ground state of 6 Li by 1.2 MeV when all 
possible angular momentum configurations are taken into account. To avoid this, we keep 
only the L — 0, S — 1 configuration. The error we make with this truncation is acceptable 
because only a ~ 5% L / component |TP| in 6 Li is neglected, but in compensation the 



asymptotic form of the wave funtion is more realistic. The contribution to I?gt from the 
6 He (L = 1)— » 6 Li transition would be negligible (~ 0.1%) since the GT operator only 
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connects states with the same L. However, because of the lack of the L = 1 component, 
the weight of the L = component is increased so that Sqt might be overestimated by 5% 
although the variational principle might partly compensate this effect. Also the presence of 
L ^ components in the 6 Li wave function might slightly shift the positions of its nodes. 
The importance of these errors can be estimated from the calculation of BQ^(g.s.). 
After these choices, our 6 Li wave functions read with L — 0, S — 1, and J n = 1 + , 

N d -1 



* 6 jm= E A{[[$ a $t} sX ??(p ad )} JM } 



i=0 



+ 



A{[[&* h ] 3 tf( Pth j\ JM } 



(7) 



for the ground state, and 



*f M = A{[[$ a ^] s gf(E,p ad )] jM } 

+ E A[[* a *i]s&( E >Pad)] JM } 

i=l 

+ A{[[^ h ] sX t L h (E,p th )] JM \ 



(8) 



for the scattering states (h = 3 He), where E is the a + d relative motion energy in the 
center of mass frame. The normalization of g is chosen consistently with as 



9o d (E,p ad ) ^Y 00 (p ad )pJ 

x (Fo(kpad) cos 5 + G {kp ad ) sin 5 



(9) 

if Pad oo, where k is the wave number, F and Go are Coulomb functions, and 5 is the 
s-wave phase shift, at energy E. 

The deuteron cluster being very distortable we use 5 basis states (N d = 5) for its de- 
scription. The ground state $q with —2.20 MeV energy and 2.1 fm (point nucleon) rms 
radius, and four pseudostates $f are included in (|7|) and Because the distortion of the 
a particle is weak in scattering processes, we describe it with a single stabilized size param- 
eter. Using the same force as for 6 He, the ground state of 6 Li described by (0) provides the 
experimental a + d cluster separation energy of 1.47 MeV. It is remarkable that we can 
reproduce both the 6 Li and 6 He ground states with the same force in this model. Notice 
that other configuration choices underestimate the binding energy, (i) Without the t + 3 He 
component, 6 Li is underbound by more than 0.5 MeV. (ii) This result is not improved by 
a distortion, (iii) With a distortion and the t + 3 He component, 6 Li is also underbound 
but by the smaller value 0.24 MeV; the role of t + 3 He is weakened because of a higher 
threshold. For the purpose of tests described below, it is interesting to have another wave 
function reproducing the experimental 6 Li energy without the t + 3 He component. This can 
be achieved by refitting the exchange mixture parameter to the value u = 0.97. 

The scattering states are calculated from a Kohn-Hulthen variational method for the 
S-matrix, which uses square integrable basis functions matched with the correct scattering 



asymptotics [[17]]. To make the calculations numerically stable, the matching radius must be 
chosen in the 10-15 fm region. In order to calculate Sgt analytically up to 25-30 fm, the 
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scattering wave function g coming from the variational method is also expanded in terms of 
square integrable tempered Gaussian functions. The squared deviation between g and this 
expansion is variationally minimized up to 40 fm and becomes less than 10 -8 . The L = 
a + d phase shifts in the E = 0.1 — 5.0 MeV relative center of mass energy region are 
compared with experiment [2(J in Fig. [I]. They agree within 0.1 degree with those obtained 
with the potential of Ref. [21], which is fitted to experiment. Putting the \I/ He , \1/ Ll , and 
ty ad (E) wave functions, coming from the variational calculations, into (Q) and we can 
compute the desired -Bqt and BQ^(g.s.) matrix elements. All calculations are performed 
analytically with the aid of a symbolic computer language. 
Comparing our model to the previous theoretical models 



1C.11.12] we can say the fol- 



lowing. A common drawback of those three models is that the a + d scattering states are 
obtained from an a + d nucleus-nucleus potential. Although this non-microscopic potential 
model reproduces the correct a + d phase shifts, the inner part of the scattering wave func- 
tions is not well established. It may not be consistent with the 6 He description. In Ref. [TI| 
6 He is represented by a pure a(nn) configuration. The dineutron component is therefore 
overestimated in the 6 He wave function, and hence in its halo part. Consequently the exter- 
nal part of the GT matrix element and Bqt are too large. In Ref. ||11|| , the results are rather 
insensitive to the use of an attractive or a repulsive interaction in the s wave, which both 
reproduce the phase shifts. This insensitivity is rather surprising because phase-equivalent 
potentials usually do not lead to the same results for off-shell effects [22|. As we can see e.g. 
in Fig. 3 of Ref. ||11|| , the repulsive potential does not produce an internal node in the a + 
d wave function, while the attractive potential produces such a node (Fig. 2 of Ref. [[[(J). 
Therefore, the GT matrix elements with the 6 He wave function should be very different from 
each other. This effect is probably hidden by the orthogonalization procedure adopted by 
the authors of Ref. [IT]. In Ref. [0], the 6 He wave function is not totally free for the varia- 
tional method: in the halo, the spatial part is restricted to small interneutron distances. As 
the interaction, employed there, does not reproduce correctly the n + n scattering and binds 
the dineutron, it might distort the halo in a larger variational space. Moreover neglecting 
the ~ 15% L ^ component in that wave function leads to an enlargement of -Bgt- Finally, 
using a microscopic wave function as initial state with an a + d potential wave function 
as final state, even with the correct phase shifts, is not consistent. The present model is 
essentially free of such defects. The description of 6 He is consistent with that of 6 Li and of a 
+ d, the interaction satisfies conditions (i) to (vi) of Sec. I and does not bind the dineutron, 
and L ^ components are explicitly included in the description of 6 He. 



III. RESULTS AND DISCUSSION 

For i?GT(g-s.) we obtain 4.60A 2 in our model which corresponds to a half life t\i2 ~ 835 
ms, to be compared with the experimental value 806.7±1.5 ms [^j. In the test model, where 
the t + 3 He component is omitted in 6 Li and the force is refitted, we get 4.48A 2 for BQ^(g.s.) 
and 857 ms for t\/2- This example shows the importance of using the same interaction in 
both the initial and final states of the decay process: the error on the theoretical result is 
doubled in the test calculation. The 3-4% deviation from experiment in the full model can 
be considered as a good agreement. Of course the lack of L ^ 1 component in 6 Li probably 
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reduces the difference between theory and experiment. Nevertheless, as discussed in Ref. ||, 
deviations as large as 7-8% could be accounted for by meson exchange current corrections. 
Such corrections are typical of different weak processes in few-nucleon systems (see e.g. Ref. 



The calculated dW/ dE curves are shown in Fig. |^ together with the experimental points 
of Ref. f9j. We can see that both the full, consistent, model (solid curve) and the test model 
(long-dashed curve) undershoot the experimental points, but the results of the consistent 
model are much closer to them. As for Bq^^.s), the deviation from experiment is roughly 
twice as large in the test model as in the full one. The error on BQ^(g.s) is typical of 
a process occurring at the normal nuclear scale, i.e. at usual distances between nucleons 
in a nucleus. This factor of two is still observed in dW/dE. However the 3-4% error on 
i?GT(g-s.) in our best model is amplified to 60-70% on dW/dE. This emphasizes again that 
the small I?gt value in the (3 decay to a + d is the result of a balanced cancellation between 
two larger values with different signs. The total branching ratio in the consistent model is 
6.0 • 10~ 6 . The branching ratio for deuterons whose energy is above the experimental cutoff 
|J is 3.1 • 10~ 6 . Let us note that the integration of the experimental data in Fig. |2| provides 
5.4 ■ 10" 6 , i.e. less than the value (7.6 ± 0.6) ■ 10~ 6 quoted in Ref. |]. 

To seek further after possible origins of errors, we changed the number of basis functions 
between 8 and 25 in the critical region of a(nn) relative motion in 6 He, each basis always 
ensuring the covering of this region up to 20-25 fm. Our results remain stable within one 
percent. Next we checked the asymptotics of our 6 He wave function. In Fig. |3] we show the 
radial part of the Xo[oo"o relative motion function at fixed p nn = 2 fm. One can see that the 
Gaussian cutoff starts only at roughly 30 fm. The correct asymptotic behavior of the bound 
state wave function of three neutral particles which have no bound binary subchannels is 
given by ]25 | 



^ -»• £- 5/2 exp(-K£) (10) 

if g-^oo. Here g is the hyperradius (defined by g 2 = Airf, where is the center-of- 
mass coordinate of cluster i with respect to the center of mass of the three-cluster system) 
and k = ^mAT-E^/fr 2 ) 1 / 2 , where Eb is the three-cluster separation energy and win the 
nucleon mass. Since we are interested in the a(nn) relative motion, g can be expressed 
as g = J f p 2 ( nn ) + \p 2 m - We checked that our Xo[oo]o re l a tive motion function satisfies the 
asymptotic form ( [T0|) (the other channels only slightly contribute when p nn is small and 
p a{nn) is large). 

Finally, to test the contribution coming from the halo part of the 6 He wave function, we 
choose to expand the p a ( nn ) space in (j^) with 10 basis functions which cover only the 0-10 fm 
region. Of course this truncation does not affect the 6 He energy at all. The corresponding 
asymptotic behavior of Xmoolo * s a ^ so sriown in Fig. [| The resulting dW/dE is the short- 
dashed curve in Fig. |^. As we can see, neglecting the 8 — 10 < p a ! nn ) < 25 fm part of the 
6 He wave function can cause an order of magnitude error. This confirms that the branching 
ratio is very sensitive to the halo part of the 6 He wave function, as was pointed out in Ref. 



12 



The theoretical curve displays a satisfactory order of magnitude although the data appear 
to be underestimated by about a factor 1.7 between 0.6 and 1.45 MeV. Its shape resembles 
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the experimental curve, but translated by about 200 keV toward lower energies. As such a 
translation does not have any obvious theoretical explanation, let us first discuss to which 
extent experiment constraints the shape of the theoretical curve. The data curve is significant 
only between 0.6 and 1.45 MeV. Below this region, the experimental cutoff hides the energy 
dependence of the results. Very few events (five) are observed above 1.45 MeV. Hence the 
real shape of dW/ dE and in particular the location of its maximum remains open questions. 
As we now discuss, these questions are of theoretical importance. 

Most theoretical results (those of the present microscopic calculation and of its two 
approximations displayed in Fig. those of Refs. pH|JTT| , and several cases studied in Ref. 
T2|| ) display a common energy behavior with a maximum between 0.3 and 0.4 fm. Only 



the magnitude of dW/ dE varies strongly. Notable exceptions are given by the EH3 and EH 
results in Ref. [Q whose maximum is at or beyond 0.5 fm, and the EH2 result which displays 
a zero near 0.6 MeV in contradiction with experiment. These non standard shapes follow 
from small modifications in the external component of the GT matrix element for a fixed 
internal component. In spite of their apparently better agreement with the available data, 
the EH3 or EH shapes cannot be considered as confirmed. Experimental data below 0.5 MeV 
(cm.) would therefore provide information about the nature of the cancellation mechanism 
between the internal and external components of the GT matrix element JT2]. They would 
discriminate between 'standard' curves with similar shapes but varying magnitudes, and 
'non standard' ones with a maximum at larger energies which is compatible with - but not 
established by - the presently available data. The standard shape is well reproduced by 
the present microscopic calculation. The 1.7 multiplicative factor might partly be due to 
the difficulty of determining the absolute normalization of experiment. The non-standard 
shapes require a very specific cancellation mechanism and would indicate the need for further 
improvements of the theory. 

The discrepancy between theory and experiment for dW j dE can possibly be attributed 
to the same mechanism as the few-percent discrepancy in i?GT(g-s.). A rather small cor- 
rection acting mainly on the internal part of the matrix element which would provide a 
correct Bq^^.s.) should modify much more strongly dW/dE and might even affect its 
energy dependence. A good candidate is provided by meson-exchange currents. Indeed 
meson-exchange corrections might have the necessary order of magnitude and should affect 
differently the internal and external parts of the GT matrix elements since the distances 
between the subsystems are different. Rather small modifications of the internal part may 
result in 50 to 100% modifications in dW/dE as discussed in Ref. ]T3], and as shown here by 



the test calculation. The fact that the tensor force is not taken into account in the descrip- 
tion of the a + d scattering may also play some role, but is probably less important than 
for the 6 Li ground state. Its introduction in the model would however exclude the use of the 
Minnesota force and require the difficult construction of a new force satisfying conditions 
(i) to (vi) of Sec. I. 

IV. CONCLUSION 

In summary, we have studied the f3 delayed deuteron emission from 6 He in a microscopic 
cluster model. We have chosen the Minnesota nucleon-nucleon interaction which reproduces 



8 



the bulk properties of the free clusters and gives very good agreement with the experimental 
N + N and a + N phase shifts in all relevant partial waves. This interaction was shown to 
provide a good overall description of the ground state of 6 He 0. In the present work we 
point out that this interaction also provides high quality results for the ground state and a 
+ d scattering states of 6 Li. All these results are obtained without any free parameter. The 
6 He, 6 Li, and a + d wave functions are used to calculate the Gamow- Teller matrix element 
of the (3 decay process. 

In a calculation of the transition to the 6 Li ground state we emphasize the importance 
of using a consistent description for the systems appearing in the initial and final states, 
respectively. The results for the (3 decay toward continuum states show the strong sensitivity 



to the halo part of the He wave function predicted in Ref. ||12|| . Our best model gives a (3 
delayed deuteron spectrum close to the experiment: the branching ratio is 3.1 • 10~ 6 to be 
compared with (7.6 ± 0.6) • 10~ 6 Q, when the experimental cutoff is assumed. The total 
theoretical branching ratio is 6.0 • 10~ 6 . We argue that the difference between the theoretical 
and experimental dW/ dE has the same origin as the small difference between the values for 
-BoT(g-S-) in spite of the fact that it is much larger. We conjecture that corrections due to 
meson-exchange currents are important enough to explain these differences. 
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FIGURES 



FIG. 1. a + d phase shift for the s wave obtained with the microscopic wave function VPj^f 
[Eq. (|8|)]. The experimental points are taken from [20|. 



FIG. 2. Transition probability dW/dE per time and energy units in the center of mass frame as 
a function of the center of mass energy E. The experimental points are taken from Ref. 0, while 



the statistical error bars are from Ref. [12|. The solid curve is the result of the present calculation 
with the t + 3 He component included in the a + d wave function, while the long-dashed curve is 
the result without this component. The short-dashed curve is obtained by dropping the p a r nr A > 10 
fm part of the 6 He wave function. 

FIG. 3. Radial part of the Xo[oo]o re l a ti ve motion function of the 6 He wave function (full line) 
calculated for p nn = 2 fm, and its truncated approximation (short-dashed curve) corresponding to 
the short-dashed curve in Fig. g. 
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